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In  agriculture,  pesticides  in  both  liquid  droplet  and  solid  particle  form  are  used  extensively,  but 
their  use  is  frequently  hampered  by  lack  of  optimum  efficiency  and  precision  of  application,  despite  the 
best  efforts  of  manufacturers  and  applicators.  This  results  in  an  excessive  economic  burden  on  agriculture 
in  the  maintenance  of  the  best  of  quality  in  agricultural  produce,  as  well  as  an  increased  ehange  of  con- 
taminating surrounding  areas.  Problems  of  air  pollution  are  also  of  prime  pertinence  and  concern.  It  is 
desired  that  the  results  of  these  investigations  in  fine  particle  behavior  and  its  inherent  subject  matter 
wiD  serve  agriculture  and  its  allied  industries,  and  other  industries  with  similar  vexations,  in  the  allevia- 
tion of  these  problems. 


2 


SOME  THEORETICAL  ASPECTS  OF  SOIL  FUMIGANT  DIFFUSION 


Ross  D.  Brazee  and  Orve  K.  Hedden  ^ 


Calculation  of  soil  fumigant  distribution  patterns  on  the  basis  of  diffusion  theory  is  con- 
sidered for  several  typical  application  configurations.  The  diffusion  formalism  of  Hemwall  is 
used,  and  techniques  for  solution  are  discussed.  Instantaneous  plane  sources  for  one  dimen- 
sion and  instantaneous  line  sources  in  two  dimensions  are  essential  to  the  formal  approach. 


INTRODUCTION 


A number  of  workers  have  studied  soil  fumigation  from  the  standpoint  of  diffusion  theory.  For 
example,  Thorne  (8),^  Siegel,  Erickson,  and  Turk  (7),  and  Lange  and  Carlson  (5)  have  made  experimental 
studies  of  fumigant  diffusion  patterns.  However,  the  papers  of  Goring  (2),  Pieczarka  and  Warren  (6)  and 
Hemwall  (3,  4)  are  of  particular  interest  to  this  discussion.  Goring  outlined  some  of  the  pertinent  dif- 
fusion concepts  and  discussed  some  experimental  studies.  The  experimental  work  of  Pieczarka  and 
Warren  dealt  with  the  effect  of  surface  sealants  on  fumigant  diffusion.  Hemwall  developed  a diffusion 
equation  with  a special  diffusion  coefficient  formulation  to  account  for  the  effects  of  soil  porosity,  tempera- 
ture, organic  matter,  water,  and  mineral  content  on  the  fumigation  process.  Hemwall  then  calculated 
dosage  patterns  for  a plane  normal  to  instantaneous  infinite  line  sources,  and  studied  the  effects  of  the 
various  soil  factors  through  their  influences  upon  the  diffusion  coefficient.  To  surmount  the  difficulties  of 
biological  dosage  computations,  Hemwall  used  digital  computer  techniques. 

Diffusion  theory  calculations  can  be  expected  to  be  useful  in  determining  where  the  fumigant  should 
be  initially  placed  in  order  to  treat  a specified  soil  region.  Hence,  the  purpose  of  the  subsequent  discus- 
sion will  be  to  study  some  mathematical  procedures  for  treatment  of  some  typical  fumigant  application 
configurations  on  the  basis  of  Hemwall’s  formalism. 

This  paper  will  consider  the  soil  layers  as  homogeneous  and  isotropic  with  a fixed  diffusion  coefficient. 
The  case  of  non-homogeneous  soils  or  layered  soils,  i.e.,  composite  media,  though  recognized,  will  not  be 
treated  at  this  time.  Extensive  use  will  be  made  of  the  unit  instantaneous  source,  either  as  an  infinite 
plane  normal  to  the  flow  direction  in  the  one-dimensional  case,  or  as  an  infinite  line  perpendicular  to  the 
flow  plane  in  two  dimensions.  These  concepts  are  illustrated  in  figures  1,  2,  and  3.  Figure  1 shows  the 
initial  configuration,  or  condition,  at  ^ = 0 later  described  {see  Sec.  7)  as  a unit  instantaneous  plane  source 
at  depth  x'  in  the  region  0 < x < a,  with  the  soil  surface  sealed  and  an  impermeable  surface  at  x = a.  A 
second  variation  is  source  and  soil  surface  alone  with  no  surface  seal  (surface  loss  of  fumigant  is  ruled  to 
occur  whenever  the  seal  is  absent)  and  without  the  impermeable  surface  at  x = a so  that  essentially  a semi- 
infinite medium  exists.  A third  variant  of  this  configuration  would  have  the  source-soil  surface  combi- 
nation with  surface  seal  absent  but  an  impermeable  barrier  at  x = a.  The  fourth  has  the  source-soil  surface 
combination  with  surface  seal  present  and  without  an  impermeable  barrier  below  at  x = a.  Each  of  the 
configurations  shown  in  figures  2 and  3 may  be  applied  in  four  similar  variations.  The  instantaneous  “slab 
source”  of  figure  2 amounts  to  a combination  of  instantaneous  plane  sources.  The  instantaneous  infinite 
line  source  is  illustrated  in  figure  3. 

The  instantaneous  source,  as  a mathematical  model,  is  intended  to  represent  the  usual  physical  mode 
of  fumigant  application.  In  practice,  a specified  quantity  of  fumigant  is  injected  at  the  desired  point  at 
t = 0,  and  allowed  to  spread.  This  is  precisely  analogous  to  the  nature  of  the  mathematical  instantaneous 
source.  It  is  only  necessary  to  multiply  the  unit  source  solution  by  the  fumigant  concentration  per  unit 
length,  area,  or  volume  as  the  case  may  be. 


’ Pioneering  Research  Laboratory  on  Physics  of  Fine  Particles  and  Pest  Control  Equipment  and  Practices  Section,  respectively. 
Agricultural  Engineering  Research  Division,  ARS,  USDA,  Wooster,  Ohio. 

^Parenthetical  numbers  in  italics  refer  to  Literature  Cited  at  end  of  publication. 
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Seal 


Soil  Surface 
(x=0) 

Sou 


Barrier 

Figure  1.  — Configuration  concept  for  an  instantaneous  plane  source  at  depth  x = x'  from  soil  surface,  with  an  impermeable  barrier 
at  depth  x = a.  Although  an  uneven  soil  surface  is  shown  in  the  drawings,  this  factor  is  neglected  mathematically;  and  in  the  case 
of  sealants,  intimate  seal-soil  contact  is  assumed. 


ce 


Figure  2.  — Configuration  for  multilocated  instantaneous  plane  sources  distributed  atx'  > 0,  0 < x'  < a,  with  an  impermeable  barrier  at 

depth  x = a. 
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Figure  3.-Configuration  for  an  array  of  instantaneous  infinite  line  sources  at  depth  x — x',  with  an  impermeable  barrier  at  depth  x = a. 


An  example  shown  in  figure  4 corresponds  to  the  condition  discussed  (^ee  Sec.  9.  I)  in  which  the  con- 
figuration of  figure  2 is  varied  by  utilizing  only  the  source  (rototilled  into  the  soil),  and  a seal  (polyethylene 
film)  placed  over  the  soil  surface.  In  figure  5,  the  method  of  injection  at  a point  behind  a traveling  blade 
(similar  to  the  “field-cultivator”  type  of  application)  corresponds  to  the  variation  of  the  figure  3 configura- 
tion {see  Sec.  12)  for  the  unsealed  soil  surface-line  source  combination  in  the  absence  of  the  barrier  at 
x = a. 


Figure  4. — Fumigation  treatment  being  rototilled  into  the  soil  and  covered  with  a polyethylene 

surface  seal. 
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The  subsequent  discussion  deals  basically  with  the  solution  of  some  boundary  value  problems,  certain 
ones  of  which  are  standard  in,  e.g.,  heat  conduction  theory.  The  surface  and  barrier  structures  will  set  the 
boundary  conditions,  and  the  apphcation  configuration,  or  instantaneous  source  arrangement,  will  estabhsh 
the  initial  conditions.  The  precise  mathematical  statement  of  the  initial  and  boundary  conditions  will  be 
covered  in  a later  section  (see  Sec.  2). 

It  is  important  to  be  aware  of  certain  advantages,  and  limitations,  of  the  diffusion  theory  formulation 
of  the  fumigation  process.  The  diffusion  equation  arises  naturally  from  random  flight  considerations 
which  could  be  expected  to  be  appropriate  to  the  process  involved.  As  such,  the  equation  represents  the 
behavior  of  some  average  property  of  a many-particle  ensemble  and,  has  a very  definite  directionality  in 
time.  From  theorems  of  thermodynamics,  irregularities  in  such  average  properties,  e.g.,  initial  instantane- 
ous sources,  will  be  smoothed  with  the  progression  of  time.  The  diffusion  equation  will  thereby  represent 
processes  where  entropy  increases  as  time  passes.  In  applying  ordinary  diffusion  theory,  one  must 
recognize  that  in  a specific  problem,  as  time  departs  from  r = 0,  the  concentration  of  the  diffusing  sub- 
stance wiU  immediately  take  on  some  finite,  though  perhaps  very  small  value  at  aU  points  in  the  space  of 
concern.  This  is,  of  course,  physically  impossible,  and  arises  from  the  fact  that  the  ordinary  diffusion 
equation  assumes  an  infinite  velocity  of  propagation.  However,  this  will  be  no  serious  difficulty  in  practi- 
cal problems,  particularly  after  a very  short  time  period. 


1.  THE  DIFFUSION  EQUATION 


HemwaU  (4)  uses  a diffusion  equation 


dc(x,y,  t) 


dt 


D 


d^c(x,y,t)  d^c(x,y,t) 
dx^  dy^ 


— krc(x,  y,  t). 


(1) 
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where  c{x,  y,  t)  is  the  concentration  of  fumigant  in  the  soil  water,  kr  is  the  rate  constant  for  irreversible 
chemical  decomposition  of  the  fumigant  in  the  soil,  and  the  unsteady-state  diffusion  coefficient  D is  given  by 

0.66kaD^ 
ka  + ku:KnKa 

In  (2),  ka  is  the  volume  fraction  of  the  soil  that  is  continuous  air  space;  Dq  is  the  diffusion  coefficient  of  the 
fumigant  through  air;  ku-  is  the  volume  fraction  of  the  soil  that  is  water;  Kh  is  the  Henry’s  law  constant 
defined  as  the  ratio  of  the  concentration  of  the  fumigant  in  the  soil  water  to  the  concentration  in  the  soil 
air;  and  is  the  adsorption  constant  defined  as  the  ratio  of  the  moles  of  fumigant  adsorbed  at  the  various 
solid  surfaces  per  unit  volume  of  soil  to  the  concentration  of  the  fumigant  in  the  air  space. 

If  the  transformation 


c(x,  y,  t)  =u{x,  y,  t)e  V, 
is  defined,  then  (1)  takes  the  standard  form 


(3) 


d~u  d'^u  _1  du 
dx^  dy^~Ddt' 

which  will  be  used  throughout  this  paper.  The  same  transformation  will  apply  to  the  boundary  conditions 
which  are  subsequently  discussed. 


2.  TYPES  OF  BOUNDARY  CONDITIONS  AND  SOLUTIONS 


The  initial  and  boundary  conditions  for  the  problems  of  concern  will  always  be  of  the  general  form 


c — Ufiixi),t^O,  a,  < Xi  < bi; 

i 

(5) 

dc 

^,-c  = 0,  x,  = Ui,  t >0,  i = 1,  2; 

uX>\ 

(6) 

dc  ' 

- — l-/3,c  = 0,  Xi  = bi,  t >0,  1,  2, 

axj 

(7) 

where  the  o;,,  al,  /Si  and  j3i  are  constants.  Hence,  the  discussion  of  Carslaw  and  Jaeger  (i,  pp.  33—35)  on 
product  solutions  of  the  diffusion  equation  may  be  expected  to  apply,  although  not  without  care.  By 
virtue  of  the  transformation  (3),  the  initial  and  boundary  conditions  will  generally  appear  in  the  form 


u =H fi{xi) , f = 0,  a,-  < Xi  < bi; 
i 

(8) 

ai^  — BiU  = 0,  Xi  — Qi,  t > 0,  i=l,2: 
axi 

(9) 

ct'i  ^3-  I3iu  — 0,  Xi  = bi,  ? > 0,  i = 1,  2, 
dxi 

(10) 

In  this  paper,  it  will  be  assumed  that  the  rate  of  fumigant  loss  through  the  soil  surface  into  the  atmos- 
phere (in  which  the  fumigant  concentration  effectively  vanishes)  is  linearly  proportional  to  the  concentra- 
tion just  inside  the  surface.  The  mathematical  statement  of  this  condition  is 

^ = hu,  x^O,  t >0,  (11) 

where  jc  = 0 is  taken  at  the  soil  surface  and  h is  a constant.  This  is  analogous  to  the  linear  “radiation” 
boundary  condition  of  the  theory  of  heat  conduction  (1,  pp.  18-19).  If  a plastic  film  or  other  sealant 
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which  may  be  considered  impenetrable  is  placed  on  the  soil  surface,  the  boundary  condition  statement 
becomes 


— = 0,  x = 0,  t>0. 
ax 


(12) 


The  vertical  downward  direction  will  usually  be  taken  along  the  positive  jc-axis  of  coordinates.  If  the 
soil  is  pervious  to  the  fumigant  over  a considerable  depth,  the  region  of  interest  will  be  x > 0.  However, 
if  the  upper  regions  of  the  soil  are  underlaid  by  a relatively  impervious  layer,  then  the  region  of  interest 
will  be  0 < ;c  < a,  and  for  the  present  discussion  we  will  presume  to  write 

1^=0,  x = a,  r > 0,  (13) 


with  the  barrier  at  depth  a. 


3.  THE  UNIT  INSTANTANEOUS  SOURCE^ 


In  one-dimension,  the  unit  instantaneous  plane  source  at  x'  normal  to  the  direction  of  flow  is  written 


2{7rDt)'^ 


(14) 


In  two-dimensional  problems,  the  unit  instantaneous  line  source  at  {x' , y')  normal  to  the  plane  of  flow  is 
written 


^ g-[{x-x')^+(y-y')^]l4Dt  ^ 

4>7rDt 


(15) 


4.  GREENES  FUNCTIONS  IN  THE  SOLUTION  OF  THE  EQUATION  OF  DIFFUSION 


In  this  section,  the  Green’s  function  method  for  solving  the  diffusion  equation  will  be  outhned.  The 
terminology  of  Carslaw  and  Jaeger  (1,  p.  356)  will  be  followed. 

Two  basic  types  of  Green’s  function  are  appropriate  to  consider  in  diffusion  theory.  First,  consider- 
ing the  two-dimensional  case,  we  may  take  the  Green’s  function,  G(x,  y,  t\x' , y' , r),  to  be  the  concentration 
of  diffusing  substance  at  (x,  y)  at  time  t due  to  a unit  instantaneous  line  source  generated  at  {x' , y')  at 
time  T,  with  the  concentration  initially  at  zero  throughout  the  region  and  the  bounding  curve  held  at  zero 
concentration.  Secondly,  the  Green’s  function  may  be  taken  as  the  concentration  at  (x,  y)  at  time  t due  to 
a unit  instantaneous  line  source  generated  at  {x' , y')  at  time  r,  with  the  boundary  normal  gradient  of  con- 

dll 

centration  specified  in  the  manner  t — = hu.  In  either  case,  it  is  found  that  the  Green’s  function  solution 

uTt} 

of  (4)  for  a region  bounded  by  a curve  s is 


w(x,  y,  t)  = j j [G{x,y,t\x  ,y' ,r)']T=Qf{x' , y')dx'dy'  + Dj 


cf){x',y’,T)  ^ds' 

oTLl 


dT, 


(16) 


where  f{x,  y)  is  the  initial  distribution  of  diffusing  material,  (f>{x,  y,  t)  is  the  surface  concentration,  and 
djdjii  signifies  differentiation  along  the  inward  normal  to  s.  For  the  one-dimensional  case,  statement  (16) 
becomes 

u(x,t)^j  [G(x,  t\x',T)]r^o/(x')dx'  + D j (f){T)^  dT.  (17) 


^See  pp.  255—259  of  Literature  Cited  1. 
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The  physical  interpretations  of  (16)  and  (17)  with  respect  to  soil  fumigant  diffusion  are  simple  but 
important:  The  concentration  at  time  t in  the  soil  region  with  initial  concentration  f{x,  y)  (or  f{x))  and  zero 
concentration  just  outside  the  surface  is  obtained  from  a distribution  of  instantaneous  sources  at  ^ = 0 
over  the  region,  with  an  amount  of  substance  f{x,  y)  dx  dy  (or  f{x)  dx)  being  liberated  in  the  element  dx  dy 
(or  dx)  at  x,  y (or  x).  Hence,  when  the  Green’s  function  for  the  given  region  and  boundary  conditions  is 
known,  the  solution  of  the  problem  for  an  initial  concentration  which  is  an  arbitrary  function  of  position 
can  be  immediately  written  on  the  basis  of  the  forms  (16)  and  (17). 

The  Green’s  function  method  will  be  found  helpful  in  several  of  the  problems  of  this  paper.  In  some 
cases,  the  problems  consist  largely  of  finding  the  Green’s  function.  We  now  proceed  to  the  treatment  of 
specific  problems,  starting  with  one-dimensional  flow. 


5.  UNIT  INSTANTANEOUS  PLANE  SOURCE  AT  DEPTH  x'  AT  t = 0 \N  THE  REGION  x>0. 

SOIL  SURFACE  SEALED 


The  appropriate  form  of  the  diffusion  equation  is 

d^u  _ 1 du 
dx^  D dt 


(18) 


with  the  boundary  condition 


|^  = 0,;c  = 0,^>0. 

dx 


The  solution  may  be  written  down  immediately  as 

1 


u(x,  t) 


|^g-(x-j')2/4Z)?_|_  g-(x+x')*/4Z)f  J 


(19) 


(20) 


2{7rDt)'^ 

by  the  method  of  images  (i,  p.  273).  By  transformation  (3),  the  solution  in  terms  of  fumigant  concentration 


IS 


c{x,  t) 


e 


2{ttDi) 


— g-(x-x')*/4rXLj_  g-(x+x')*/4I)f  J , 


(21) 


6.  UNIT  INSTANTANEOUS  PLANE  SOURCE  AT  DEPTH  x'  AT  1=0  IN  THE  REGION  x>Q. 

SOU  SURFACE  UNSEALED 


The  boundary  condition  is 


^ = /iM,  JC  = 0,  ^ > 0. 

dx 


(22) 


Carslaw  and  Jaeger  (7,  358-359)  outline  the  Laplace  transformation  solution  in  detail.  The  result  is 


u{x,  t) 


1 


2i7TDty 


^ g-(x-x')*/4ni_j_  g-(x+x')2/4Z)l  j j^^Dth‘^+h(x+x')Qj.j^Q 


2{Dty 


(23) 


which  for  fumigant  concentration  becomes 


c{x,  t)  = e->'d  I + ^-(x+xnv^oq  _ f^^Dth2+h(x+x’)erfc 


^ X I 


(24) 
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In  (23)  and  (24),  the  expression 


erfc 


is  the  error  function  complement. 


_ 2{Dty 


w 


+ hVUt 


(25) 


7.  UNIT  INSTANTANEOUS  PLANE  SOURCE  AT  DEPTH  x'  AT  ^ = 0 IN  THE  REGION  0<x 
< a.  SOIL  SURFACE  SEALED  AND  AN  IMPERMEABLE  SURFACE  AT  x = a 


The  boundary  conditions  are 


We  first  write 


where 


du 

dx 


= 0,  x = 0 and  x = a,  t > 0. 


U = V-\~W, 


1 


2{7TDty 


(26) 


(27) 


(28) 


and  ic  is  a solution  of  the  one-dimensional  diffusion  equation  which  approaches  zero  as  t approaches  zero 
and  is  chosen  so  that  the  boundary  conditions  (26)  on  u are  satisfied. 


The  Laplace  transform  of  (27)  is 
and  for  (28)  is 

where 

The  subsidiary  equation  for  W is 


Z{u}  = U=V-\-W, 

y—  -I.  ..  p-q\x-x'\ 

2qD 

q-ipIDf. 


(29) 

(30) 

(31) 


dW 

dx^ 


qW=0. 


Hence,  U is  given  by 


U — - sinh  qrjc  + 5 cosh  qjc, 

2qU 


(32) 


where  A and  B are  arbitrary  constants.  Upon  applying  the  boundary  conditions  (26)  and  simplifying  the 
results,  we  find  that  U takes  the  form 


and 


rr  cosh  q{a-\-x  — x')-\- cosh  q{a  — x~x')  „ ^ ^ , 

U = ry; r-r , 0 < .x  < AC  , 

Zuq  Sinn  qa 


cosh  qf(a  + x' —x)  + cosh  g(a  — x' —x)  , 

U— ,x  <x<a. 

Zuq  sinh  qa 


(33) 


(34) 
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The  inversion  theorem  may  be  apphed  to  determine  u{x,  i),  viz., 


u{x,t)=S^  '^{U{x,p)]  = 


y+iR 

hm  I 1 r e^Tcosh  q{a  + x — x')  +cosh  q{a  — x — x')^dp 

1 27ri  J 2Dq  sinh  qa 

y — \R 


(35) 


where  (33)  is  to  be  considered  first.  The  integrand  in  (35)  has  simple  poles  at 

p = Q;  p = — DnVla^,n=\,2, . . . (36) 

We  select  a contour  of  the  form  shown  in  figure  6 and  large  enough  to  contain  aU  the  poles  of  the  integrand. 
The  residues  at  the  poles  of  the  integrand  are  (i/a)  for  p = 0,  and 


n f e^^[cosh  q{g-\-x  — x')-\-  cosh.q{a  — x — x')^  1 

1 2Dqsmhqa 


= innla 


Q-DrChrHla^ 


2D[{dldp)  iipID)’^  sinh  (a(p/i))^))] 


r.  2 2,  2 

p=-Dn  7T  la 


e {a-\-x  — x )n7T  , {a  — x — x )mr 

; COS h cos  

a cos  {mr)  a a 


(37) 


1] 


Then  the  sum  of  the  residues  inside  the  contour  will  be 


S Res  = - + - y 

« «n=l 


TITTX 


COS 


mrx 


(38) 


We  have  for  the  integral  on  the  contour  of  figure  6 

1 \ f ^~^”e^^[cosh  g(a  + jc:  — :r')  + cosh  q{g  — x — x'y]dp 
2m  1 J y+ioo  2Z>g  sinh  qa 


= y Res, 


(39) 


since  the  integral  along  the  large  circle  is  zero.  Hence,  for  the  inverse  transform  we  have 


u{x,  t)  = - 
a 


1 + 22 

n=l 


e-Dn^^tla^COS 


( mrx 
V a 


cos 


mrx 


(40) 


An  identical  result  is  found  for  the  transform  (34),  which  indicates  (40)  is  a complete  solution  in  terms  of 
u{x,  t). 

The  solution  to  our  problem  in  the  notation  of  fumigant  concentration  is,  therefore. 


c{x,  t) 


71=1 


COS 


/ mrx  \ / mrx'  \ 


(41) 


8.  UNIT  INSTANTANEOUS  PLANE  SOURCE  AT  DEPTH  x'  AT  1 = 0 IN  THE  REGION 
0<:t;<  a.  SOIL  SURFACE  UNSEALED  AND  IMPERMEABLE  SURFACE  AT  x = a 


The  boundary  conditions  are 


— = hu,x=Q^t>Q; 

(42) 

^ = 0,  X = o.  ^ > 0. 
ox 

(43) 

The  solution  may  be  obtained  by  Laplace  transform  methods,  but  we  shcill  merely  adapt  the  general 
result  of  Carslaw  and  Jaeger  (7,  p.  360)  and  write 


u{x,t)  = '^  Zn{x)Zn{x')e-^^¥, 
n=l 


(44) 


where 


Zn{x) 


(2)^o:n[Q:n  cos  janx)  + h sin  {anx)  ] 
{{a^  + h?)  {aal  + han\ Y 


(45) 
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and  ±a„,  «=  1,  2,  . . . , are  roots  of 

aiaxi  {aa)  =h.  (46) 

Carslaw  and  Jaeger  (7,  p.  491)  give  tables  of  the  first  six  roots  of  the  transcendental  equation  (46). 

The  solution  in  terms  of  fumigant  concentration  is 

c{x,  t)  =e-V^Zn{x)Zn{x')e-^^nt.  (47) 

7*=1 


9.  SOLUTIONS  FOR  ARBITRARY  INITIAL  DISTRIBUTIONS  OF  CONCENTRATION 


In  obtaining  the  solutions  for  u{x,  t),  in  Sections  5 to  8,  we  have  at  the  same  time  obtained  Green’s 
functions  G{x,  t\x\  r)  for  the  regions  of  concern.  Hence,  it  is  possible  to  write  down  a series  of  solutions 
for  arbitrary  initial  distributions  of  fumigant  by  using  (17)  and  the  appropriate  Green’s  function.  Mathe- 
matically, the  initial  condition  wiU  be 


c{x,  0)  = f{x) , 0 < ^ < a or  00 
The  solutions  for  several  regions  foUow. 


(48) 


I.  The  region  x>  0 with  the  soil  surface  sealed.  Boundary  condition  of  Section  5. 

y-krt  C » 


c{x,  t) 


e 


2{TTDtf 


(49) 


II.  The  region  x>  0 with  the  soil  surface  unsealed.  Boundary  condition  of  Section  6. 


c{x,  t)  = e-'^vj  le-^x-x')mDt  + ^-(x+x')mDq  _ f^^Dth2+hix^x')erfc 


x-\~x' 


2{Dty 


^h{Dt) 


f{x')dx'.  (50) 


III.  The  region  0<x<  a with  the  surface  sealed  and  an  impermeable  surface  at  x=a.  Boundary  con- 
ditions of  Section  7. 


p-k  tf°  r » 
<;(;£,«)  = — 1 + 22. 
^ 7o  L 


eos  ff)  COS 


f{x')dx' . 


(51) 


IV.  The  region  0 < x < a with 
conditions  of  Section  8. 


the  surface  unsealed  and  an  impermeable  surface  at  x=a.  Boundary 


c{x,t)  = e I f Zn(x)Zn(x')f(x')dx' 

n~T  w . 


(52) 


The  expressions  for  Zn(x)  and  otn  are  as  given  in  (45)  and  (46),  respectively.  The  vahdity  of  the  interchange 
of  the  operations  of  integration  and  summation  has  been  assumed. 
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10.  TWO-DIMENSIONAL  CONFIGURATIONS 


(15), 


W e consider  a unit  instantaneous  line  source  atx  = x'  and  y = y/ . 


This  source  may  be  written  as  in 


— i — p-l(x-x'P+(y-y':)^]/4Dt 

47TDt  ^ 


(53) 


With  the  help  of  this  general  form  and  the  Green’s  functions  derived  in  the  previous  eases,  we  will  attempt 
to  write  down  some  solutions  for  the  two-dimensional  line-source  array.  It  is  assumed  in  all  cases  that  at 
great  distances,  the  concentration  approaches  zero.  Positive  or  negative  integers  are  assigned  to  j as 
y'j  is  positive  or  negative,  and  in  the  order 

y-m<  . . . <y-^  <y-i  < 0 <yT<y2'<  . . . y'n.  (54) 


11.  UNIT  INSTANTANEOUS  SOURCE  ARRAY  AT  DEPTH  x'  AT  t=0  IN  THE  REGION  x>0, 

SOIL  SURFACE  SEALED 


From  the  boundary  condition 


1^  = 0,  x = 0,  i > 0, 

dx 


and  equations  (20)  and  (53),  we  write  down  the  solution 

C{X,  y,  t)  [e-^x-x'mDt^  g-(x+x')2l4Dq  e-(y-y^2i4Dt 

4-7701.  . 


j = -M 


(55) 


(56) 


12.  UNIT  INSTANTANEOUS  SOURCE  ARRAY  AT  DEPTH  jc'  AT  ^ = 0 IN  THE  REGION  x>0, 

-a.<y<(X).  SOIL  SURFACE  UNSEALED 


By  the  boundary  condition  (22),  and  equations  (23)  and  (53)  we  write 

c{x,  y,  t) 


^ ^ I 1 ^g-(x-x')^l4Dt  _|_  g-(x+x')2/4Dr  j 


2{7rDtr.^_^  [2(7tZ)0 


' h{Dtf 


_2iDtf 


e-iy-yp^l4Dt.  (57) 


13.  UNIT  INSTANTANEOUS  SOURCE  ARRAY  AT  DEPTH  AT  ^ = 0 IN  THE  REGION 
0<x<a,  ^<y<o).  SOIL  SURFACE  SEALED  AND  AN  IMPERMEABLE  SURFACE  AT 

x = a. 


From  the  boundary  conditions  (26),  and  equations  (40)  and  (53)  we  have 

-kj 


c{x,y,  t) 


e ’'r 


2a{7rDt)'^ 


S j 1 + 2 2;  cos  1 — 1 cos 


Q-(y-y^)‘^l4Dt  ^ 


(58) 


14 


14.  UNIT  INSTANTANEOUS  SOURCE  ARRAY  AT  DEPTH  a;'  AT  t = 0 IN  THE  REGION 
0<x<a,  ^ <y  < 00.  SOIL  SURFACE  UNSEALED  AND  AN  IMPERMEABLE  SURFACE 

M x = a 


On  the  basis  of  the  boundary  conditions  (42)  and  (43),  and  equations  (44)  and  (53)  we  obtain 

c(x,  y,  0 = 2^^  i;  [ 5 Zn(x)Zn{x')e-‘>-^'^  (59) 

with  Zn{x)  and  a,i^  as  expressed  by  (45)  and  (46),  respectively. 
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